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Abstract 

Estimation of leaf productivity in eelgrass {Zostera marina L.) is crucial for evaluating 
the ecological role of this important seagrass species. Although leaf marking 
techniques are widely used to obtain estimates of leaf productivity, the accuracy of 
these assessments, has been questioned mainly because these fail to account for leaf 
growth bellow the reference mark and also because they apparently disregard the 
contribution of mature leaf tissues to the growth rate of leaves. On the other hand, 
the plastochrone method is a simpler technique that has been considered to 
effectively capture growth in a more realistic way, thereby providing more accurate 
assessments of both above- and below-ground productivities. But since the actual 
values of eelgrass growth rates are difficult to obtain, the worth of the plastochrone 
method has been largely vindicated because it produces assessments that 
overestimate productivity as compared to estimates obtained by leaf marking. 
Additionally, whenever eelgrass leaf biomass can be allometrically scaled in terms of 
matching leaf length in a consistent way, the associated leaf growth rates can be 
also projected allometrically. In this contribution, we used that approach to derive an 
authentication of the plastochrone method and formally demonstrate that, as has 
been claimed to occur for leaf marking approaches, the plastochrone method itself 
underestimates actual values of eelgrass leaf growth rates. We also show that this 
unavoidable bias is mainly due to the inadequacy of single-leaf biomass assessments 
in providing a proxy for the growth of all leaf tissue in a shoot over a given interval. 
Moreover, the derived formulae give conditions under which assessments of leaf 
growth rates using the plastochrone method would systematically underestimate 
matching values obtained by leaf marking procedures. And, assessments of leaf 
growth rates obtained by using the present data show that plastochrone method 
estimations underestimated corresponding proxies obtained allometrically (27%), or 
through leaf marking (35%). Allometric projection is recommended as a simpler and 
more effective procedure to reduce the bias in eelgrass leaf productivity estimations 
that associates to the use of plastochrone methods. 

Keywords: Eelgrass, Allometric scaling, Leaf growth, Plastochrone method, Formal 
validation 



BioMed Central 



© 2013 Echavarria-Heras et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Echavarria-Heras et al. Theoretical Biology and Medical Modelling 2013, 10:34 
http://www.tbiomed.eom/content/1 0/1 /34 



Page 2 of 15 



Introduction 

Eelgrass (Zostera marina) is a cosmopolitan seagrass species that plays an important 
role in shallow and nearshore ecosystems. This temperate macrophyte is distributed in 
Northern Hemisphere habitats from the Arctic Circle to the Tropic of Cancer [1], 
where it provides a nursery for fish and as a substrate for attached algae and epi- 
fauna [2,3]. By fixing large amounts of carbon through photosynthesis, eelgrass 
plays an important trophic role, sustaining detrital food chains and other secondary 
producers [4], Eelgrass also helps in the remediation of contaminated sediments 
[5] by filtering and retaining nutrients from the water column [6] and contributing 
to the stabilization of sediments [7]. Moreover, eelgrass meadows reduce erosional 
forces by stumping wave energy, thus promoting the stabilization of adjacent 
shorelines [8,9]. 

The variability in eelgrass biomass constitutes a dynamic link between its structural 
and trophic roles, because changes in the amount of organic carbon that can be fixed 
modulate the structure of the habitat for the associated biota. These organisms are af- 
fected in different ways when changes in biomass occur seasonally or unpredictably 
[10]. Therefore, accurate measurements of the standing crop and productivity of eel- 
grass constitute an important input for evaluating the ecological functions and values 
of this significant seagrass species [11]. 

Several methods have been developed to measure seagrass growth and productivity 
[12-14]. Procedures for measuring productivity in aquatic macrophytes have included 
the assessment of changes in biomass over a growing season, as well as the quantifica- 
tion of oxygen production or incorporation of 14 C during photosynthesis. But the ac- 
curacy of these methodologies was questioned early on [15-18]. Growth in seagrasses 
occurs through the expansion of modules formed by rhizome segments, which have 
bundles of attached leaves and roots. Because every leaf produced corresponds to the 
production of a rhizome node, it is reasonable to assume that seagrass growth and 
leaf formation are equivalent processes [19]. This conspicuous feature has encour- 
aged efforts to estimate the growth of eelgrass, as well as that of other seagrasses 
with ribbon-like leaves, by measuring leaf growth. Leaf growth in seagrasses can be 
estimated by using the leaf marking technique which was originally proposed as a 
way of avoiding the underestimation of growth when using oxygen metabolism 
measurements [16,18]. Leaves were marked above the sheath with a small staple at 
the height of a reference frame placed above the sediment, and then recovered 
after a period of 2-4 weeks. The new growth in each leaf between the reference 
frame and the staple was then weighed. The ratio of this weight to the number of 
days of growth determined the leaf growth rate. Because this procedure is consid- 
ered to give consistent estimations [19,20] it has been used and modified by sev- 
eral authors (e.g. [21-23]). 

For assessment of growth in Zoostera marina in particular, workers have marked 
blades using felt- tip pens, staples and hypodermic needles [4,24-26], with the reference 
level being either the leaf tip [25,27] or a point at the top of the sheath [12,24,26,28]. 
However, Brouns [29] asserted that this technique only estimates blade production and 
is not suitable for measuring total production. The reliability and accuracy of the leaf- 
marking method has been also questioned by others; for example, because it does not 
account for new leaf growth within the sheath below the reference mark [24], or it 
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might inflict damage to tissues [19], or fail to account for leaf maturation processes 
[11,13,30], or that it requires destructive dry- weight measurements. 

Based on the observation that the production of every seagrass leaf is associated with 
the production of a rhizome node, Patriquin [31] envisioned that seagrass root and rhi- 
zome production can be also estimated by determining the time interval for the forma- 
tion of new leaves and by counting the associated leaf scars. The time elapsed between 
the appearance of two consecutive leaves has been termed a plastochrone interval or 
plastochrone index [32,33]. For eelgrass, Jacobs [4] used the average weight of the third 
leaf on sampled shoots as a substitute for the leaf biomass gained by a shoot over a 
plastochrone interval and divided this leaf biomass surrogate by the value of the 
plastochrone interval to assess growth rates of leaves. Gaeckle and Short [11] endorsed 
the method of Jacobs [4] and used the weight of a mature leaf to represent all growing 
leaf tissue in a shoot over a given a plastochrone interval to calculate eelgrass growth, 
calling the leaf-biomass to plastochrone ratio the plastochrone method for eelgrass 
leaf-growth assessments. And Gaeckle and Short [11] concluded that leaf marking re- 
sults in lower estimates of leaf growth than the plastochrone-based method, as asserted 
by Brouns [29], and they considered that this approach both fully captures growth and 
is also simpler and non-destructive. 

But despite the advantages of plastochrone methods, their use has not been yet for- 
mally substantiated. In the present research we attempt to fill this gap by using allomet- 
ric models and a discrete mathematical formulation for the increment in biomass 
gained by an eelgrass leaf over a time interval to derive an equation which expresses 
the mean shoot leaf-growth rates in terms of a factor of the leaf biomass to 
plastochrone ratio plus a remainder, and explore their analytical implications. And, we 
extend aforementioned result to provide a formal device aimed at obtaining the abso- 
lute deviations between leaf growth-rate assessments estimated by using the conven- 
tional leaf-marking technique and those obtained by using the plasthochrone method. 
We also corroborate the derived formulae using both simulation and real data. Finally, 
we stress the advantages of using allometric methods in eelgrass research and discuss 
the findings of this study. 

Data and related calculations 

The data used for this study were collected in a Z marina meadow in Punta Banda Estu- 
ary, located near Ensenada, Baja California, Mexico. Form March 1999 to July 2000 we 
visited the site biweekly. At each sampling time using the Kentula and Mclntire technique 
[26] approximately 40 shoots were marked, and those previously marked were retrieved. 
On all the shoots collected, the sheaths were peeled off and the leaf contents were sepa- 
rated. We measured the lengths of each of the retrieved leaves and the associated leaf- 
length increments gained over the marking interval. The matching biomasses were 
obtained by direct dry weight determinations. The associated in situ values for the mean- 
shoot leaf growth over the marking interval were estimated directly by dividing total leaf 
production in the marked shoots retrieved by the number of days elapsed. 

Matching estimations, produced by the leaf marking method, were obtained by 
subtracting leaf production within the sheath to the in situ rates. Total leaf productivity 
in retrieved shoots was also determined indirectly by the allometric approach described 
in equation (12a) below, and also by applying the plastochrone method. 
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A formal derivation of the plastochrone method 

In this section we show how, using allometric models for the representation of eelgrass 
leaf biomass in terms of length, we can provide a formal substantiation of the 
plasthocrone method for leaf-growth assessments. With that aim, we use the index 5 to 
label a generic Zoster a marina shoot and let n(s) stand for the number of leaves it 
holds. This number includes wholly formed leaves attached to the sheath and those de- 
veloping inside this structure. We let lj s (t) stand for the length of j-th leaf in the 
addressed shoot, and use At to denote a positive amount of time. Thus, the symbol Alj s 
(t, At) stands for the increment in length attained by the leaf lj s (t) along the interval 
[t, t + At], i.e. 

Al js {t,At) = l js (t + At)-l js (t) (1) 

Correspondingly, the symbol o)j s (t) with 1 <j < n(s) denotes the dry weight at a time t 
of the j-th leaf in the considered shoot. Also, we will let Aodj S (t, At) denote the gain in 
dry weight for l JS (t) occurring along the interval [t, t + At], i.e., 

Aco Js (t, At) = co Js (t + At)-(o js (t) (2) 

Lets denote by means of the symbol Lg s (t, At) the mean leaf growth rate attained by 
a shoot s during the interval [t,t + At]. We then have the mean shoot leaf-growth rates 

s-^AcoJt.At) , . 

Lg s (t,At) = Y, ' (3) 

Similarly the symbol Lg(t, At) represents the average of the Lg s (t, At) values taken 
over the number of shoots collected at a time t + At, i.e., 

L ^ = ^TW (4) 

Where N(t + A£) is the number of shoots collected at a time t + At. In what follows 
the calculated values for the Lg(t, At) rates will be referred to as in situ mean leaf- 
growth rates. 

Assuming that there are parameters a and such that coj s (t) and lj s (t) are linked 
through the allometric equation 

o js (t) = al js (tf (5) 

then a related expression for the leaf biomass increment A(dj S (t, At) can be derived [34] 
and this becomes 

Ao) ajs (t,At) = al js {t + Atf 8j 8 (t,At) (6) 
with 

And, since we have 

0<Al Js (t : At)<l Js (t + At) (8) 
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we must also have, 

0<5 Js (t,At)<l (9) 
Since the result of equation (6) is derived from equation (5) we must formally have 

Aco js (t,At) =A(o ajs (t,At) (10) 

Nevertheless, for applications to a given data set, when calculating, the leaf biomass 
increments of equation (6) we must rely on estimates for the parameters a and hence 
the associated uncertainties will set, 

AcD js (t,At)=Aco aJs (t,At)+Rj s a (t,At) (10a) 

where RJ s a (t,At), is a residual term. 

Correspondingly we will denote by means of the symbol Lga s (t,At) the allometric sur- 
rogate of the Lg s (t, At) rates of equation (3). That is, 

L g a s (t,A t) =Y p^ (11) 

Again we must formally have, 

Lg s (t,At)=Lga s (t,At) (11a) 

and, as it was elaborated in equation (10) when estimates of the allometric parameters 
are used to calculate the Ao) a j S (t,At) proxies we must have, 

Lg s (t,At)=Lga s (t,At)+R™(t,At), (lib) 

where R™ a (t,At) is the involved estimation error. 

Averaging the Lga(t, At) values over the number of shoots collected at time t + At, 
one gets allometrically projected values Lag(t, At) for the in situ leaf growth rates Lg(t, At) 
of equation (4), these are 

Lga s (t.At) 

W,*) = % + At) • (12a) 

Furthermore we must formally have, 
Lg(t,At) =Lga(t,At) +Ba(t,At) (12b) 

where the term Ba(t, At) stand for the pertinent approximation bias. Again, only in the 
case in which both equations (10) and (11a) hold we will have that the Lg(t, At) rates of 
equation (4) and their allometric representations Lga(t, At) are equivalent, this setting Ba 
(t, At) = 0. In what follows the result of equation (12a) will be called "allometric method 
for the projection of eelgrass leaf growth rates" or simply as "allometric method". 

Now let p(t, At) denote the value of the leaf plastochrone interval calculated using the 
leaf marking data resultant for the interval [t,t + At] [4], and lets assume that the biomass 
of the third leaf on each retrieved shoot provides a substitute for all growing leaf tissue in 
that shoot along a time interval of size p(t, At) [4,11]. Then, in accordance with our con- 
vention, the biomass of this leaf will be denoted by means of G) 3s {t,At), and lets also rep- 
resent the resulting leaf biomass to plastochrone ratio, using the symbol Lgp s (t, At) that is, 

^-ot (i3) 
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this equation is the formal representation of the plastochrone method, which is consid- 
ered to yield the mean rate of growth for all the leaves on a retrieved shoot s and over the 
interval [t, t+At] [4,11,14]. Then, averaging the shoot rate values Lgp s (t, At) over the 
number of shoots collected at a time t + At one obtains the proxy values Lgp(t, At) pro- 
duced by the plastochrone method and corresponding to the in situ leaf growth rate 
values Lg(t, At), they are, 

L ^ At) = N(t + At) ■ (14a) 
Again by letting Bp(t, At) stand for the associated bias we have, 

Lg(t, At) = Lgp(t, At) + Bp(t, At) (14b) 
Now, assuming that equation (10) holds then from equations (11) and (13) one obtains 

Lg s (t,At) = C? (t,At)l4p,{t,At)+S?(t,At), (15) 
where the factor C ap (t,At) is given by, 

(?M,|=«f*, (16) 
and R a s p (t, At) is a positive remainder defined through, 

jy (Mt) = ^^^ + ^W' 3 ) (17) 
with, 



|l otherwise 



Since G)j S {t, At) is positive and we have also observed that 8j S (t, At) does not simultan- 
eously vanish for all values of the index ; in a retrieved shoot, we must have, 

R ap (t,At) >0. (19) 
Also, from inequality (9) and equation (16) we obtain, 

0<Cf (t,At)<^^- (20) 

Then, since depending on the order relationship that p(t, At) and At satisfy, the ratio 
of p{t, At) to At can take any positive value, and since as given by equation (16) Cf 
(t,At) depends continuously on S 3s (t } At) } p(t, At) and At it can formally take any posi- 
tive value. For instance we could set C ap (t,At) = 1, which according with equation 
(16) leads to the equivalent equation 

s ^ At) =wk) (21) 

which taking into account inequality (9) could be only satisfied whenever the order re- 
lationship At < p{t y At) holds. 
Moreover, from equations (4), (13) and (15) one obtains, 

Lg(t, At) = Cap(t, At)Lgp(t, At) + Bap{t, At) (22) 
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where 



Cap(t,At) 



J2 s Lgp s (t,At) 



(23) 



and 



Bap{t,At) 



(24) 



N(t + At) 



so, inequality (19) implies, 

Bap(t,At) > 0, 



(25) 



and using inequality (20) and equation (23) we get, 



0<Ca P (t,At)< P -^. 



(26) 



Again, since C a s p (t,At) varies continuously and Lgp s (t, At) does not vanish and de- 
pends continuously on both to 35 (£, At) and p{t, At) from equation (23) we conclude that 
Cap(t, At) varies continuously and takes all the positive values within the range set by 
inequality (26). 

From equation (22) we can formally infer that the equivalence of Lg(t, At) and Lgp 
(t, At) would occur whenever the constrains Bap{t, At) = 0 and Cap(t, At) = 1 simultan- 
eously hold. Moreover from equation (23) we have that Cap(t, At) = 1 is equivalent to 
the equation, 



which will be satisfied whenever C a s p (t,At) = 1 and which according to the elaboration 
sustaining equation (21) can only occur whenever the order relationship At < p(t, At) 
holds. Indeed from inequality (20) we have, 



then if we assumed that At >p(t,At) the factor ( C ap (t,At)-l) in equation (27) would 
not change signs, and since Lgp s (t, At) is positive for all values of t equation (27) could 
not be satisfied. But even though, as we have elaborated above, whenever At < p{t, At), 
the statement Cap(t,At) = l Cap{t, At) = 0 could be formally satisfied, inequality (25) 
implies that the equation Bap{t, At) = 0 could not be satisfied, precluding the equiva- 
lence of Lg{t y At) and Lgp(t, At). As a result, when projecting the Lg(t, At) rates by 
means of the Lgp(t, At) proxy values, we can anticipate a systematically positive or 
negative bias. This direct bias has been denoted by means of Bp(t, At), (cf. eq. 14b) and 
from equation (22) it is given by 




(27) 



-l<Cf (t,At)-l< 



p(t,At)-At 
AtT~ 



(28) 



Bp[t,At) = [Cap[t,At)-\)Lgp[t,At) +Bap{t,At) 



(29) 



Then, since as given above Bp(t, At) approaches the positive value Bap(t, At) whenever 
Cap(t, At) approaches one, and since Bp(t, At) varies continuously and does not changes 
signs, it must remain positive in all its variation range, then as a result Lg(t, At) values 
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will be systematically underestimated by their Lgp(t, At) surrogates. Moreover, using in- 
equality (26) we can find bounds for the variation of Bp(t, At) that is, we have 

BLap(t,At)<Bp(t,At)<BLap(t,At) + P ^^ Lgp(t,At) (30) 
where 

BLap(t,At) = Bap(t,At)-Lgp(t,At) (31) 

This result shows that the minimum possible value that Bp(t, At) can attain is BLap 
(t, At) and that whenever the inequality p(t, At) > At holds the upper bound for its vari- 
ation range will increase from the value Bap(t,At). Moreover this order relationship will 
induce wider variation ranges for Bp(t, At) relative to the complementary case p(t, At) 
< At, therefore substantiating the claim by Short and Duarte [14] who stated that for an 
efficient application of the plastochrone method the length of the observation period 
At must exceed the value p(t, At) of the leaf plastochrone interval. 



Derivation of a formal relationship between leaf marking and plastochrone 
method assessments of eelgrass leaf growth rates 

We now show that the addressed allometric framework can also be used to obtain a 
formal connection between leaf marking technique assessments and those obtained by 
means of the plastochrone method. With that aim lets now denote by means of the 
symbol Lgrn s (t, At), the estimations of the Lg s (t, At) rates of equation (3) which are 
obtained by means of the leaf marking technique. Lets also denote by means of Lgm(t, At) 
the average of the Lgm s (t, At) values over the number of shoots collected at time t + At 
that is, 

^ Lgm s (t,At) 

Lgm(t, At) = ^ s . . (32a) 

N(t + At) K J 

And, if Bm(t, At) denotes the related approximation bias we formally have 

Lg(t,At) =Lgm(t 1 At) + Bm(t,At). (32b) 

Now, defining the ratio 

_ Al js (t,At) 

x ' s ^ m) -WTm) — (33) 

then, since Al is positive and bounded above by lj s (t + At),we have, 

0<X Js (t,At)<l (34) 
and, from equations (6) and (33) we get, 

Aco ajs {t,At) =a(Alj s (t,At)Y( 1 ^- Xlit ' At) A (35) 

V Ai\t,Aty ) 

which rearranging leads to 

Ao ajs (t,At) = p(Al js (t,At)) a +R™(\,t,At) (36) 
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where 

Rf s a {\, t, At) = pi Js (t + Atf (l-(l-A/(f , At)f-X js (t, Atf^j (37) 

Hence, Ad) a j S (t,At) can be split into two components: a term (>(Alj S (t, At)) a ; the alio- 
metric expression of the biomass of the leaf length increment A/, plus a remainder, R™ a 
(A, t,At). The term ft(Alj S (t, At)) a , can be linked to the contribution of new leaf material 
to growth. And, since R^ a (X,t,At) } depends on the whole leaf length span lj s (t + At), 
this residual term can be taken as the contribution of mature leaf material to growth. 
Moreover, formally we must have 

ijm<Mt)=E ««r. (38) 

/ 

and assuming that equation (10) holds, equations (3) and (36) imply 

Lg s (t,At) = Lgm s (t,At)+R™(t,At), (39) 
where 

RT a (t, At) = J2 ' At (40) 

therefore, equations (15) and (39) yield, 

Lgm 8 (t,At)= Cf (t,At)Lgp s (t,At)+R™P(t,At) (41) 

where 

R^(t,At) = -£ ^ t (42) 

j 

then, averaging over the number of shoots collected at time t + At we have that the 
Lgrn(t, At) and Lgp(t, At) rates are related through, 

Lgm(t,At) = Cap(t,At)Lgp(t,At) + Bmp (t,At), (44) 
where 

-S^lr <«) 

then as it was the case in equation (22), whenever the inequality, 

Bmp(t,At)>0 (46) 

holds, the result of equation (44) implies that the Lgrn(t, At) assessments will be 
underestimated by the Lgp(t, At) values. 

Corroboration of derived results using the present data 

Using the present leaf data, we calculated the values for in situ leaf growth rates Lg(t, At) 
by means of equation (4). The concomitant Lgrn(t, At) estimations produced by the 
leaf-marking technique were also obtained. Using the formula in Jacobs [4], we pro- 
duced biweekly estimations of plastochrone-interval values. Plastochrone method 
proxies Lgp(t, At) for the in situ leaf growth rates Lg(t, At) were also calculated. The 
model of equation (5) was identified as consistent with the present data set, producing 
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a coefficient of determination of 0.82. The values of the fitted parameters were 
a = 0.00002 and ft = 1.3 similar to those previously reported for this site [34,35]. 
We then used these parameter values and associated leaf data to produce allometric estima- 
tions Lga(t, At) for the in situ leaf-growth rates using equation (12a). An analysis of variance 
found significant differences between mean annual leaf-growth rates among methods 
(F= 17.95 df: 3, \\%p< 0.001) and an a posteriori Tukey test found non-significant differ- 
ences among in situ values and those produced by the leaf marking method (p = 0.89), 
among in situ values and those obtained through the allometric method (p = 0.08), and 
among leaf marking values and those produced by using the allometric method p = 0.29. 
The Tukey test also showed that the mean annual leaf-growth rate obtained through the 
plastochrone method is statistically different from that linked to the in situ values, that calcu- 
lated from leaf-marking data, and the one obtained from allometric projections (p < 0.001). 

Figure la provides a direct comparison of measured Lg(t, At) leaf growth rates and 
those calculated, by using the Lgrn(t, At), Lga(t, At) or Lgp(t, At) methods. It is shown 
that all these proxies produced positive bias in approximating the values for in situ Lg 
(t, At) rates, but Bm(t, At) was smaller relative to Ba(t, At) and at the same time this bias 
was found to be smaller than Bp(t, At). Moreover, since the bias term Bmp{t, At) was also 
found to be positive, then equation (44) implies that the values of the plastochrone- 
method proxies Lgp(t, At) underestimate the Lgm(t, At) values obtained by means of the 
leaf-marking technique. Differences on the root mean squared error (RMSE) values on 
Table 1 suggest that while the leaf marking and allometric proxies could be expected to 
produce consistent projections of observed eelgrass leaf growth rates, the plastochrone 
method instead shows a relatively smaller reproducibility of the named rates. 
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Table 1 Root mean squared error (RMSE) and concordance correlation coefficient (p) 
values for comparison of observed leaf growth rates and corresponding proxies 
calculated by the leaf- marking, allometric or plastochrone methods 



Criteria 


Leaf marking 


Allometric 


Plastochrone 


RMSE 


0.00021283 


0.000586175 


0.001700075 


P 


0.9895 


0.9831 


0.5286 


95% Confidence interval 


(0.9776, 0.9951) 


(0.9638, 0.9921) 


(0.2010, 0.7498) 



Moreover, an analysis of Lgrn(t, At) versus Lg(t, At) values found a statistically significant 
regression through the origin; Lg(t,At) = 1.040 Lgrn(t, At), r 2 = 0.98,/? < 0.001 (Figure lb), 
also when we analyzed Lga(t, At) versus Lg(t, At) the same regression model was found 
highly consistent; Lg(UAt) = 1.163 Lga(t,At), r 2 = 0.96, p < 0.001 (Figure lc). Meanwhile, 
when we compared Lgp(t, At) versus Lg(t, At) no regression through the origin could be 
fitted. Instead only a poor fit; Lg(t, At) = 0.779 Lgp(t, At) + 0.002, ? = 0.26, p < 0.003 could 
be produced (Figure Id). Additionally, values of the Lin [36] Concordance Correlation 
Coefficient (p) for reproducibility of the Lg(t, At) rates by the methods addressed here 
(Table 1), reveal that both the leaf marking (p = 0.9895) and allometric (p = 0.9831) 
methods performed better than the plastochrone method (p = 0.5286). Indeed we can 
observe that the bias Ba(t, At), produced by the allometric method of equation (12a) 
is relatively smaller than that linked to plastocrone method assessments. Since the 
allometric approach is a simpler technique than leaf marking procedures [34,35], 
allometric projection is recommended as an effective method to reduce the bias as- 
sociated to plastochrone method assessments. 

Discussion 

Allometric scaling models of the form Y - aX^, where Y is a response and X is an 
explanatory variable, are commonly used in biological research. In this formulation, a 
is positive and is known as the normalization constant, while is called the scaling 
exponent. From an empirical standpoint these allometic models have been fitted to 
many pairs of biologically traits X and Y in a highly consistent way, providing reliable 
methods for estimating variables that are difficult to quantify by using other variables 
that can be easily and directly measured. Empirical scaling models can be viewed as 
static relationships in which both the allometric exponent and the normalization constant 
take values that pertain to a particular data set. But beyond a purely empirical advantage, 
allometric approaches have demonstrated that prediction is not an unachievable aim in 
ecology. Indeed, allometric methods can be successfully used to predict roles and charac- 
teristics of organisms [37,38]. The predictive realm of allometric laws, moreover, is not 
circumscribed only to individual properties of organisms but can also apply to energy 
flows in whole ecosystems [39]. For example, the consistency of allometric relationships 
between body size and metabolism, which span 27 orders of magnitude in body mass for 
living organisms in both aquatic and terrestrial ecosystems, has provided a basis for the 
development of the metabolic theory of ecology [40] . 

The benefits of allometric methodologies in seagrass research were stressed by 
Duarte [41], who exemplified the use of static allometric models to grasp the implica- 
tions of differences in plant size for productivity. For eelgrass in particular, Hamburg 
and Homann [28] and Echavarria-Heras et al [35] used empirical allometric methods 
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for the assessment of aboveground biomass. But allometric methodologies in eelgrass 
research are not limited only to empirical approaches for nondestructive appraisals of 
leaf biomass, but have also been used to formally substantiate empirical paradigms 
employed in eelgrass research [42]. They demonstrated that allometric scaling can pro- 
vide criteria for evaluating the suitability of the leaf-biomass-to length ratio as a para- 
digm for nondestructive leaf biomass evaluations. And from a theoretical standpoint, it 
has also been established that allometric models can be linked to a paradigm for the in- 
direct assessment of leaf growth rates in terms of simple measurements of leaf length 
or area [34]. 

It has been considered that conventional leaf-marking techniques produce biased es- 
timations of eelgrass leaf-growth rates [11]. One fundamental objection states that, 
since these methods measure leaf growth by accounting only for the biomass of all the 
pieces of leaf material that have developed during the period between the marking and 
recovery of shoots, any leaf material produced below the reference point is not cap- 
tured, and thus that portion of leaf production is missed. But even though the marking 
point could be placed on the sheath itself, as a way to avoid this inconvenience, the 
modified method still relies on newly produced leaf tissue, which is considered to be 
immature and flaccid. These tissues have a lower weight-to-length ratio than mature 
leaves [11], resulting in an underestimation of leaf growth. 

The plastochrone method for eelgrass leaf growth rate estimations [4,11] was pro- 
posed in order to avoid the lack of representativeness of the leaf marking techniques. 
Unlike the leaf-marking techniques, the plastochrone method has been considered to 
produce simple, direct, nondestructive and unbiased estimations of eelgrass growth. 
Nonetheless, the reliability of this method has been largely vindicated at an empirical 
level, and practitioners have justified its application by asserting that this paradigm fully 
captures the growth pattern of eelgrass, in which new plant parts are incrementally 
produced while young tissue is simultaneously maturing. 

The implementation of the plastochrone method is centered on the idea that the bio- 
mass of a single leaf representing mature tissue on a shoot, divided by the leaf 
plastochrone interval, provides a reliable substitute for the contribution of the shoot to 
the mean growth rate of leaves. But so far, a mechanistic explanation is still missing as 
to how a single leaf can capture the contribution to growth produced by all the leaves 
in a shoot over a given observation period. This is where allometric approaches can 
contribute to the formal substantiation of the method. Indeed, equation (15) was 
obtained by starting from the basic balance equation (3), then representing the contri- 
butions to growth of all the individual leaves on a shoot using the allometric represen- 
tation of equation (6); and finally, by averaging over the number of retrieved shoots, we 
obtained equation (22). This equation relates the observed values of eelgrass leaf- 
growth rate Lg(t, At) in a direct dynamical way to their proxy plastochrone method 
values Lgp(t, At), thereby providing a formal authentication for this paradigm. 

The balance equation (3) shows, moreover, that any fixed value of Lg s (t, At) links all 
individual contributions to leaf growth in an essential way, with the biomass gained by 
a particular leaf being lost by the remaining ones. Hence, a single leaf would capture all 
the leaf biomass gained by a shoot over a given growing period, if and only if the bio- 
mass of the residual leaves in the shoot remained static, which is unsupported by the 
observed growth mode of eelgrass. This is indeed what inequality (19) properly 
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expresses, and its validity formally explains that we can expect a downhill bias when 
plastochrone method calculations are used as proxies for observed leaf growth rates. 
This was in fact corroborated with the present data; the values obtained for the derived 
formulae confirmed that the plastochrone-method proxies underestimated the leaf- 
growth rates observed. 

Likewise, the gains obtained from the use of allometric approaches are also shown by 
equation (37), which provides a device for the estimation of the contribution of leaf 
maturation processes to growth. But using the present theoretical settings, we can also 
address the controversy between the leaf marking and plastochrone methods for leaf- 
growth assessments. Indeed, equation (44) relates leaf-growth rate estimations obtained 
by means of the leaf marking technique to those obtained through the plastochrone 
method. And, the derived formulae set conditions under which it could be expected 
that plastochrone method approximations of eelgrass leaf-growth rates would systemat- 
ically underestimate matching values obtained by using leaf-marking procedures. 
Indeed the condition set by inequality (46) was verified for the present data; then as a 
consequence, the plastochrone method proxies underestimated similar values obtained 
by means of leaf-marking methods. That is the inequality Bp(t, At) > Brn(t, At) holds 
through time, as shown in Figure 1. Further, using the present allometric framework, 
we formalized the claim of Short and Duarte [11] who stated that, for an efficient appli- 
cation of the plastochrone method, the length of the observation period At must exceed 
the leaf plastochrone interval p(t, At). And, although this condition was satisfied for the 
present data, we obtained Ba(t, At) < Bp(t, At) when approximating observed values. 
That is, the bias produced by the allometric method of equation (12a) was found to be 
smaller than that linked to plastocrone-method estimates. Furthermore, for this data 
the order relationship Bm{t, At) < Ba{t> At) < Bp(t, At) was maintained through time. It 
can be also observed from Figure 1 and from the high value of the associated Lin [36] 
Concordance Correlation Coefficient (p) for reproducibility of the Lg(t, At) rates 
through Lgrn(t, At) values (Table 1), that the contributions of leaves within the sheath, 
and of leaf maturation processes to growth, were not significant enough to reduce the 
effectiveness of the leaf marking technique. 

And although a larger bias is associated with allometric projections than with leaf- 
marking techniques, the allometric approach is recommended over the plastochrone 
method because it is simpler and more precise (Table 1). Moreover our results support 
the conjecture that, taking into account the contribution to growth of all the leaves on 
an eelgrass shoot, as is done in both the leaf-marking technique and the allometric 
method provides a better proxy than the weight of a single leaf, as traditionally done in 
plastochrone method assessments. To this, we add that the lack of representativeness 
of this single-leaf- weight surrogate could explain the deviations shown in Figure 1. 

Most methods used to assess leaf productivity in eelgrass have been vindicated at an 
empirical level. In particular, the plastochrone method has been favored over leaf- 
marking techniques because it is supposed to eliminate the alleged downhill bias linked 
to these assessments. This study confirms that a formal interpretation of assessment 
procedures is a crucial step in appraising their reliability, and we have shown how allo- 
metric methods can pave the way to achieve this important step. Formally, the 
plastochrone method will unavoidably produce biased estimations of the actual eelgrass 
leaf growth rate values. Empirical corroboration validated this inference, endorsing the 
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view that allometric scaling relationships are not circumscribed only to empirical de 
scriptions of the linkage between variables. Rather, they can be used to produce theor 
etical tools aimed towards the clarification of relevant issues in eelgrass research. 
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